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Explicit density expansions of non-equilibrium probability distribution functions for molecular 
Brownian particle in ideal gas are obtained in original form what visually implies (is exact solution 
to) the previously established dynamical virial relations. Role of these relations in unbiased analysis 
of molecular chaos properties in many-particle statistical mechanics, including the mobility 1/f noise, 
is newly investigated in clear terms of conditional probabilities and averages. 



PACS numbers: 05.20.Dd, 05.20.Jj, 05.40.Fb, 05.40.Jc 



I. INTRODUCTION 

Here, we continue consideration of special but prin- 
cipally important problem of non-equilibrium statisti- 
cal mechanics: factual statistics of random walk of 
molecular-size Brownian particle (BP) interacting with 
atoms of thermodynamically equilibrium ideal gas [H-[3| ■ 
Namely, using a formal trick like suggested in our re- 
cent work 0, we shall derive explicit series expansions 
of the BP's and many-particle probability distribution 
functions over mean density of the gas. We shall see that 
these expansions, being represented in terms of the ir- 
reducible correlation (cumulant) functions, give unique 
formal solution to the exact "virial relations" [IKS 11 0] 
(if considered as equations for the correlation functions) . 
Then we shall discuss some important consequences of 
these results, first of all, concerning existence of 1/f -type 
BP's mobility fluctuations and essentially non-Gaussian 
long-range statistics of BP's path (dramatic failure of the 
"law of large numbers" in true many-particle statistical 
mechanics). 



II. DENSITY EXPANSION OF PROBABILITY 
DISTRIBUTION FUNCTIONS 

Let N ^> 1 ideal gas atoms are contained, along with 
the BP, in a box with volume f2 = N/n > a 3 N where 
n and a denote mean density of the gas (concentration 
of atoms) and characteristic radius of BP-atom interac- 
tion, respectively. The interaction potential, $(/o) , is 
assumed repulsive and short-range. Let at initial time 
moment, t — , full normalized [N + 1) -particle distri- 
bution function (DF) of this system is 



N 



D N (t = 0) = F™(R,V)H 



E{r 3 -R) G m {vj) 
J n E(r-R)d^r 



(1) 



Here R , V and rj , Vj are coordinates and velocities of 
BP and atoms, Fjj n = F (t = 0) is BP's initial proba- 
bility distribution (normalized to unit), 

G m [v) = (2nT/m)- 3 / 2 exp(-mv 2 /2T) is Maxwell 
velocity distribution for particle with mass m , and 



E{p)=exp[-*(p)/T\. 

Clearly, if F* n (R,V) = W m {R)G M {V) (with M be- 
ing BP's mass) then the distribution ([T]), D^it — 0) oc 
W m (R) exp (-ifW/T) (with ffW being full system's 
Hamiltonian) , represents thermodynamically equilibrium 
system, but it is not statistically equilibrium (stationary) 
because of its spatial inhomogeneity in respect to BP's 
position if described by more or less localized initial dis- 
tribution W m (R) , e.g. W m (R) = 6(R) (which will be 
taken in mind below). 

At t > 0, D N {t) = exp(L( N H)D N {t = 0), where 
denotes full system's Liouville operator (corre- 
sponding to H( N ' and taking into account interactions 
with the box walls). For our system we can write 



L 



N 



(2) 



where operator Lq represents BP itself, while operators 
Lj with j > represents j -th atom itself plus its mutual 
interaction with BP. 

Further, introduce sequence of (s + l)-particle ( s = 
0,1,2 ... N) marginal probability distribution functions 
(DF) for BP and s atoms: 



F s (t) = n s f ... f d 

Js + 1 JN 



N(t) 



(3) 



= n s 



s+l 



N 



e L(Nh D N (t = 0), 



where Jj symbolizes integration over variables of j -th 
atom. Thus, in respect to atoms' coordinates, these DF 
are "normalized to volume" , i.e. defined standardly, like 
in the mentioned references and originally by Bogolyubov 
in Q. Consider them under the thermodynamic limit, - 
when iV — > oo and D, — > oo (boundaries of the container 
disappear at infinity) under fixed N/ ft = n = const , - as 
functions of the mean gas density n . 

With this purpose, it is convenient, firstly, to intro- 
duce, like in [5(, "coherent product" o of operator ex- 
ponentials by definition as follows: 



.+A S 
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Then, secondly, to use identities 

N N 

e L (JV) t = e L t Y[o[l + {e L ^ - 1)] = e Lat + e Lot o[e^'i* - 1 ] . . . o [e L ^* - 1 ] (4) 

i=l k=l l< 31 <...<j k <N 

and 



e i (JV) t _ e (L +L 1 + ...+L s )t 



[] o[l+(e^* -1)]= (5) 

j = S+l 

N-s 

= e (L +L 1+ ...+L B )t {l+ J- o[e L ^ -l}...o[e L ^ -1}} 

k=l s+l<j 1 <...< Jk <N 

Substituting this expansions to ([3]) we see that any finite-number term of corresponding series expansions of F s (t) 
is well defined in the thermodynamic limit. Combining EqsfTJ [3] and [SI after quite obvious algebra, one obtains 

F (t,R,V) = S(t)Ft(R,V) , (6) 

where the BP's "propagator" S(i) is presented by series 

00 k k r 

S(t) = +^^^11 h[e L ^ -l]g{ Xj ) = (7) 
fe=i K - 3=1 J i 

= e Lot + nJ^[e {Lo+Ll)t -e Lot }g{ Xl ) + 

+ JJyL + L 1+ L 2 )t _ e (Lo+ Ll )t _ e (Lo + L 2 )t + e L t ] g ^ g ^ + ^ 

with xj ee {pj^j} , pj ee Tj - i?, g(x) ee exp [-$(p)/T] G m (v) = E(p) G m (v) , and J. . . . = J d 3 ^ / d 3 vj 

At that, after change of variables from rj to pj = rj — R and transition to the thermodynamic limit, the components 
of full Liouville operator ^ can be chosen as 

Lq = -VV - fV P , (8) 
Li = (V — vj) V, + <J>'( Pi ) (V PJ - V P ) , (9) 

where = d$(p)/dp, V = <9/<9i?, Vj = 9/3pj , V P = d/dP, V P] = d/d Pj , P = MV and ft- = mv 3 are 

momenta of BP and atoms, and we introduced external force / acting onto BP. 

Notice that fc-th term of Q represents all such BP's trajectories what include > fc collisions with fc atoms, at 
least one collision with any of them. Indeed, if BP was not interacting with j -th of them, i.e. Lj was replaced by 
(V — Vj)Vj , then integration J\ would turn into zero any whole fc-th term of with k > j . 

Similarly, from EqsfTJ [3] and [5J it follows that 



F s (t) = e 



L +L 1 + ...+L,)t 



!oo s+k f, | s 

1 +EtT IT /°[e^-l]5(%) Fttlgfa) 
fe=i ' j=s+i J i ) j=i 



oo k 1 s+/c „ I 
e (L +L 1+ ...+L,)t pin +e (L +L 1 + ...+L s )tJ2 7 ±- ) J] / o[e^ 4 - 1] I F™ k , (10) 

fe=l ' I j=s+l •''J ' 



pin = F in -Q s(a . .j 



Formulae ©-© and (|10[) give explicit formal series expansions of non-equilibrium (time-dependent) DF in respect 
to mean gas density. 
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III. CUMULANT DISTRIBUTION FUNCTIONS AND DYNAMICAL VIRIAL RELATIONS 

In addition to the DF F s (t) = F s (t, R, V, x\ . . . x s ) , let us consider functions C s (t) — C s (t, R, V,x\ . . . x s ) defined 
by Cq = F and, at s > , the series as follow: 



C s {t) = e 



3=1 y k=l j=s+l i 



k 




00 k r 

I'-' l, eh0t U ^ -^^J) K n (11) 



(at k = , of course, there is no integration). 
From this definition it is clear that, first, C s (t = 0) = at s > . 
Second, 

F 1 (x 1 ) = F <?(£!) + Ci(zi) , (12) 
F 2 (xi,x 2 ) = F g(xi) g(x!) + Ci(xx) g(x 2 ) + g(xi) Ci(x 2 ) + C 2 (xi,x 2 ) , 

and so on, where for simplicity only atoms' arguments are exposed. Hence, C s represent irreducible (s + 1) -order 
correlations between BP and s atoms and can be named "cumulant functions" (CF). 

Third, considering them as functions of the mean gas density too, C s — C s (t, R, V, x\ . . . x s ; n) (and correspond- 
ingly F s = F s (t, R, V, x\ . . .x a ; n) ), we can write 



dC s 
dn 



C s+1 (13) 



s+l 



That are the "virial relations" which for the first time were found in an original way in [(| 0] for BP in (generally 
non-ideal and dense) fluids and then, in t 9] and PHI]) derived directly from the BBGKY equations (in some of 
these references the CFs were designated as V s ). In essence, more reasonable name for such kind of results may be 
"dynamical virial relations" (DVR) [5| . 



IV. CONDITIONAL PROBABILITIES AND 
AVERAGES, AND DVR IN THEIR TERMS 

The DVR establish important connections between 
density dependence (and hence time and space depen- 
dence) of the BP's probability distribution, F , and de- 
gree and character of BP-gas statistical correlations. The 
first of relations (IT3l) already was carefully discussed in 
0, 0], 0, [IH3], but here we shall refresh it and make 
more transparent by exploiting the useful concepts of 
conditional probabilities and conditional averages. 

Namely, let us introduce conditional two-particle DF: 

Fi{t,R,V,x;n) _ n(t,x\R,V;n) 
F (t,R,V;n) n ' 1 } 

where n(t, x \R, V; n) = n(t, p, v \R, V; n) is conditional 
mean value of instant local density (concentration) of 
atoms in the p -space under given position and velocity 
of BP. That is 

n(t,x\R,V;n) = {n{t,x)) R y <n (15) 

with 

n(t,x) 5(x -xj(t)) = (16) 

3 



= Y^5{p-{r 3 {t)-R{t)))5{v- V] {t)) 

3 

and rj(t), Vj(t), R(t), V(t) being instant values of the 
system's variables. The condition is reflected by the 
subscript of angle brackets which in turn symbolize, as 
usually, averaging over given statistical ensemble (deter- 
mined by F s m ) . 

In terms of conditional statistical characteristics the 
first of DVR (fT3|) , after its multiplication by n , clearly 
takes form 

d In F (t,R,V;n) 
d In n 

= J ' [n(t,x\R,V;n) — ng{x)]dx = (17) 

= J [n(t,p\R : V;n)~nE(p)}d 3 p , 

where we introduced also instant local gas density in con- 
figurational space, 

n{t,p) = [n(t,x)d 3 v=Y / S(p~Pj(t)), (18) 

J 3 

and n(t, p\R, V, n) = (n(t, p)}R,v,n as hs conditional 
ensemble average value. 
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If we are interested mainly in spatial dependence of 
F = C , i.e. in the BP's "diffusion law" 

W(t,R;n) = J C {t,R 7 V;n)d 3 V , 

then, evidently, instead of (fT7|) we can write 

9 In W(t,R;n) 



n(t, p \R; n) 



5 In n 

nE(p) ] d a p = 
= AN(t\R,n) 



(19) 



with n(t, p\R;n) = ( n(t, p) being conditional mean 
value of instant gas density at distance p from BP under 
given BP's position, and AN(t\R, n) corresponding con- 
ditional mean change of total number of atoms in BP's 
surroundings. 

Notice that the factor nE(p) is nothing but value of 
the n(t, p\R\n) under complete statistical equilibrium, 
when there are no statistical correlations between atoms 
and BP except those coming from their direct dynamic 
interaction. Therefore, integral in (Till)) , AN(t\R,n) , 
characterizes specifically non-equilibrium "excess" ( "his- 
torical" BP-gas correlations. 



V. DIFFUSION LAW AND FIGHT OF 
INTUITIONS 

Let us consider the BP's diffusion law, W(t, R; n) , as- 
suming that (i) the BP-atom interaction potential is re- 
pulsive and short-range, (ii) the external force is absent, 
/ = , i.e. our system all time stays in equilibrium, and 
(iii) (>t, where r is BP's mean free-flight time (or 
velocity relaxation time), r ~ (A/Vb)yl + M/m , with 
A being BP's mean free path and and Vo = y/T/M 
characteristic thermal velocity. According to standard 
reasonings of gas kinetics and probability theory, firstly, 
A = (tto ri) with a denoting characteristic radius of 
BP-atom interaction. Secondly, W(t, R; n) has Gaussian 
long-time asymptotic, 



W(t,R;n) -> W^\t,R) = 



exp(-R 2 /4Dt) 
{AirDtf/ 2 



(20) 



V 2 r 



where D = D(n) is BP's diffusivity, D(n) 
VoXy/l + M/m cx 1/n for not too dense gas (at 
47ra 3 n/3<l). 

This asymptotic expresses the celebrated Bernoulli's 
"law of large numbers" and is an integral part if com- 
mon intuition cultivated by probability-theoretical way 
of thinking in statistical mechanics. However, it never 
was deduced from rigorous statistical mechanics in itself, 
without "art of conjectures" (see comments and refer- 
ences in SiSSEMl). Now, let us combine it with 
rigorous relation (|19[) and with the inverse proportional- 
ity D(n) oc 1/n. This yields 



dlnW^K 3 r2 



(21) 



which, in turn, implies that AN(t\R, n) can take ar- 
bitrary large negative values. In other words, even in 
spite of equilibrium nature of the BP's diffusion, statisti- 
cal correlations between gas and BP can involve arbitrary 
many atoms. 

This conclusion, however, is in contradiction to another 
common intuitive notion. Namely, that of Boltzmannian 
"molecular chaos" which usually serves as logical reason 
of the "law of large numbers" . Thus, the DVR visually 
reveal falsity of both these notions. 

To avoid the contradiction, we should refuse hypothe- 
sis that rigorous statistical mechanics always allows real- 
ization of "molecular chaos" and the "law of large num- 
bers" in transport phenomena. And, hence, we should 
accept profound conjecture that BP-gas correlations have 
substantial but bounded value. That is the quantity 
AN(t\R,n) , - which can be named e.g. "correlation 
number" , - satisfies 



AA_ < AN(t\R,n) < AN^ 



(22) 



where AN- and A_/V + are positive quantities indepen- 
dent (under long-range asymptotic) on R and t and, 
expectedly (at a 3 n < 1), also on n . 

This means @, H, 0] that instead of (|20|) we have to 
write, asymptotically, 



W(t,R;n) -> 



m{R 2 /ADt) 



E(R 2 /V H 2 ) (23) 



with D = D{n) cx 1/n , where, because of (|19p and 
function &(z) possesses power-law tail: 

«4( z>1 )' a = | + AN - > 



(24) 



and 3(0) = 1 . The cut-off function S(-) (fast decreas- 
ing at infinity) can be omitted except when calculating 
higher-order statistical moments of R. Since at that, 
evidently, AN + = max z d In W/d\n n =3/2, the ex- 
ponent a can be expressed as a = AN + + AN- . Ex- 
plicit example of such (realistic) "diffusion law", - for 
probe atom of slightly non-ideal gas in the role of BP, - 
with the characteristic "correlation number" A7V_ = 2 , 
is presented in flij . Thus, probabilities of large devia- 
tions of BP's path from its typical values ( \R\ ~ y/6Dt) 
are giant as compared with predictions of the (idealistic) 
Gaussian law. 

In [1, H, 0| : these consequences of the first of DVR (fT3|) 
were considered in different fashion, in terms of charac- 
teristic volume, fl c , occupied by BP-gas correlations in 
the p -space. It appears in the above formulae if we notice 
10] that, due to non-negativity of n(t, p \R, n) , equality 
(|T9| implies inequality 



9 In W(t,R;ri) 
9 In n 



> 



n Q r 



(25) 
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(we neglect contributions ~ a 3 in comparison with Sl c ) . 
This coincides with (j2"2")l at fl c = AN_/n. An unam- 
biguous formal definition of fl c was done in @, @ ; 0|- 

The asymptotic (|2U1) - (|24I) means that BP has no a pri- 
ori certain diffusivity. Instead, a variety of possible BP's 
trajectories is characterized by some effective distribu- 
tion, U(D ;t,n) , of its diffusivity, so that 



W(t,R;n) 



(t, R) U{D' ; t, n) dD' (26) 



According to (j2"3"j) and (pM)) , the distribution must have 
definite power-law tail and thus look e.g. like 



U(D;t,n) (x 



D 



AN- 



£)AJV_+1 



-1M-^W£) (27) 



D 



V H 



with Do = Dq(ji) oc 1/n being characteristic diffusivity 
scale and H(-) nearly the same cut-off function as in 
Currespondingly, in 



r(a) 



f 



r(a-3/2) (l + z) Q 



(28) 



with a = AA_ + 3'2. Coefficient T(a)/T(a - 3/2) 
is normalizing factor ensuring that (asymptotically) 
/ W(t, R; n) d 3 R = 1 . This is direct analogue of dis- 
tributions found in [l4j (for probe atom of weakly non- 
ideal gas) and in [l5| (for BP different from non-ideal gas 
atoms [la]) in the "collisional approximation" of BBGKY 
equations 

According to the universal "generalized fluctuation- 
dissipation relations" (see e.g. and [HI, EH, E3| and 
references therein), the representation (|21)|) - (|2"T1) . with 



W { D °\t,R) 



exp [-(R- Df /Tf/ADt] 
(4ttD<) 3 / 2 



(29) 



in place of (|20|) . is valid also for non-equilibrium Brown- 
ian motion under non-zero external force, / ^ , at least 
for sufficiently small one, /A <C T . 

The resulting distribution predicts and describes 1/f- 
type fluctuations of BP's mobility fi = D/T , as well as 
at f = similar 1/f - typ e fluctuations of BP's diffusivity 
0, |, H, 0, E3, E3-E3, E3411 , thus confirming logics of the 
Krylov's "art of conjectures" [25| . 

At D Q (f/T) 2 t > I (when drift contribution to BP's 
path R(t) exceeds diffusive contribution) this mobil- 
ity/diffusivity 1/f -noise visually manifests itself in prob- 
ability distribution of R(t) 's projection onto / 's direc- 
tion @: according to (j2"o| and (|33]l. its shape nearly re- 
produces that of /7(Z3; i, n) (analogous distributions of 
charge carriers' mobilities in electronics were detected by 
time-of- flight measurements [221] V 



VI. HIGHER-ORDER CONDITIONAL 
PROBABILITIES AND VIRIAL RELATIONS 

It is not hard to see that expansions (fTP)) and (fTTj) or, 
equivalently, formulae (TT2l together with the DVR (fT3")) 



imply relations 
dn 



/ [F s+ i -g(a; s+ i)F s 



(30) 



After dividing by _F S and multiplying by n this yields 
many-particle analogue of ([T7|> . 



=2EV a (t|J2,V; a: W ) n) = 

am n 

[n(£,x|i?, V> (s) ) -g(:r)n]ete 



(31) 



with = {xi...a; s } and n(t, x \R, V, x^ , n) being 
conditional average of n(t, x) under condition that BP is 
in given state and, besides, some s atoms already occupy 
s exactly known states. 

Then, introducing many-particle DF in configurational 
space, 

W s (t,R,p^;n) = J F s d 3 Vd 3 Vl ...d 3 v s , 
where = {p\ . . . p s } , one can transform (|3"01) - (|3"Tj) to 

= AN s (t\R,p^,n) = (32) 

am n 

= J [n(t,p\R,p^,n)-E(p)n] d 3 p 

with n(t, p\R, p( s ',n) — (n(t, p)) Rp ( 3 ) n representing 
average under condition that s atoms are known already 
to be at given distances from BP. 

To interpret these formulae correctly, recall that, - by 
statistical definition of DFs F s (see also [1| and below), 
- different arguments xj of F s , relate to different atoms. 
Therefore, s atoms entering into the conditions in fact do 
not enter into n(t 7 x) and n(t, p) . As the consequence, 
a contribution they give to AN(t\R, n) in (TT91 is absent 
in AN s (t\R, p( s \ n) in (f32j) . so that we can write 



AN s (t\R,p (s \n) > AN(t\R,n) 



(33) 



At that, of course, AN s (t\R,p("\n) -> AN(t\R,n) 
when points p^ are far from the "correlation volume" 
f2 c , i.e. p^ — > oo , and equality in (|33l) can be achieved 
only at small emough belonging to fl c . 

Together with these reasonings, formulae (|3"2"j) . (j3"3")l 
and (|22|) naturally do prompt that for s > , - in addi- 
tion to (|22|) . (f23|) and (|28p . - the following asymptotical 
estimates take place (at / = ) : 



min AN s (t\R,p^,n 



-AN_ 



(34) 



W 8 (t,R;n) -> , (35) 

, , T(a s ) $(z) 

a s = AA_ + s + 3/2 , (37) 

where the minimum is taken over R and , and 
V^ s (i,i?;n) means result of suitable smoothing of 
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W s (t, R, p( s ) ;n) over small enough p^ from the "corre- 
lation volume". At that, J W s (t, R; n) d 3 R = 1 , as the 
above basic definition of W s does require. 

Thus, the sequence of functions W s (t, R;n) qualita- 
tively coincides with sequence W s (t, R)/(l + sm/M) 
from [l5| (see also [6|), and coincidence would be also 
quantitative if the "correlation number" AN— was equal 
to 1 + Mj m (see below) . This observation confirms that 
exact DVR ([30)) - ([32]) and the "collisional approximation" 
[l^ - fTBj are rather close approaches. 

One can notice that the sequence W s = W s (t, R;n) 
satisfies (at Dal jn ) simple recurrent relations 



W s+1 = 



1 



AiV_ 



_ — 

s dn 



(38) 



(more generally, nd/dn may be reolaced by —Dd/dD 



Iterating them, after some combinatoric algebra one can 
obtain 



s-l 



w a = u 



AN. 



s + n O On 



d/dr 



k=0 



AN. 



E 

k=0 



s r(A7v_) d k 



k J t(An_ + k) n dn k w ° 



(39) 



Right-hand side here is mere result of "normal order- 
ing" of multiplication and differentiation operators n 
and d/dn from left side. 

At / ^ (in the linear drift response regime), corre- 
spondingly, 



/■oo 

W s -> / W^\t,R-D'ft/T)Us(D';t,n)dD' 
Jo 

Us(D;t,n) = r{AN _ - - DAN „ +1+8 exp \- — j B 



where effective diffusivity distributions U s obey recurrent relations 



U s+ i(D;t,n) 



AN. 



d_ 

dn 



U s (D;t,n) 



at Dq oc 1/n . 



(40) 
(41) 

(42) 
(43) 



VII. GENERATING FUNCTIONALS AND VIRIAL RELATIONS 

For many purposes it is convenient to accumulate all DFs in a single generating functional (GF), 

F{t,R,V,1>\ n} = F + f] K /■•■/ F > #Bi)...^(iB.) (44) 

s=l S - Jl Js 

with ip( x ) being formally arbitrary probe function. Following the original general Bogolyubov's construction of such 
functionals in 8J (where probe function u(x) — nipix) was used instead of ip(x)), we can introduce our one as 
thermodynamic limit of 



N 



F N {t,R,V,t,n) = F +J2 



s\{n - s y.n s j x 



F s ip(xi) . ..ip(x 3 ) 



[...[ D N (t)l[[l + 1>(xj)] = [■■■[ cxp|^ln[l 

Ji Jn j=1 Ji Jn \ j=1 



i Jn 



+ 4>{x 3 )] ) D \ ( / ! = 
exp I I n(f, x) In [1 + ip(x)] dx }■ D N (t) 



where f2 = N/n, and we involved the "microscopic gas density" n(t, x) — ^2j = iS(x ~ Xj(t)) whose statistics is 
determined by the total DF .Djv(i) in turn determined by (JXJ) and ©. The latter expressions highlight statistical 
meaning of of the functional (|44l) , T = lim Tn , allowing us to write [J] : 

F{t, R, V, tp; n} = ( S(R- R(t)) S(V - V(t)) exp { / n(t,x) In [1 + ip(x)] dx 
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= F (t,R,V;n) \ ex Pjy ri(t,x) ln[l +tp(x)]dx j J (45) 

Here, as above, (. . ) Ry. n is symbol of conditional averaging under given BP's variables, so that the latter angle 
brackets represent conditional characteristic functional of random field n(t, x) , with = In [1 + ip(x)] in the role 
of probe function. 

Next, introduce, in full analogy with (01} , GF of the CFs (TITj) : 

C s ^{ Xl )...^(x s ) = P{t,R,V, (l+tp)n} , (46) 



C{t,R,V,i>;n} = C Q +^ ^ /.../ 

s =i s ' Jl Js 

V{t,R,V, X } = £ 1 \ I ■■ I x{xi).-.x(xs)e Lot f[o[e 

s=0 S - [ Jl Js j=l 



l]g{ Xj ) } Ft(R,V) (47) 



The latter equality in (|46l) directly follows from CF's definitions (ITlT) . with functional V{t, R, V, x} representing 
unified "propagator" for all CF's. In view of (fl2j) . 



F{t, R, V, ip ; n } = C{t, R, V, i/j ; n } exp \ n / i/j(x) g(x) dx } (48) 



Since, again by the CFs definition, C s (t = 0) = at s > and therefore C{t — 0, R, V, ip', n} = Fq u (R : V) , one 
can see that the exponential in (|48|) is conditional characteristic functional of thermodynamically equilibrium gas: 

1 { ' = "' '"' ] - exp/n / i>(x) g(x) dx } = ( exp ( / n(t,x) In [1 + i>(x)] dx \ ) * (49) 



F™{R,V) r U n ' sw J \ J/i^n 

Factual independence of this functional on {R, V} just reflects the gas' equilibrium, which is indicated by superscript 
"eq". The change of variable ip(x) = exp [£(&)] — 1 transforms (|4T?)) into 



e X p{n/^)[e^-l]^} = (e Xp {/n(M)^)^ 



(50) 

R,V,n 



which shows that corresponding statistics of random disposition of gas atoms in p -space is trivial Poissonian. 

However, the equilibrium of gas itself does not mean that of the whole system, if F™{R, V) differs from £1~ 1 Gm{V) . 
Then at t > the above mentioned excess "historical" correlations between gas and BP will appear to be accumulated 
by GF (|46p . For the first time these correlations, in their connections with BP's diffusion law, were considered in 
|26l I27J . In our present case of BP in ideal gas, these connections, expressed by DVRs (fT3|) . together form generating 
DVIljll, [3, H (quite obvious from g7])) 

C{t, R,V, a + ip ; n } = C{t, i?, V, ip/(l + a) ; (1 + a) n } , (51) 

where a = const is independent on x . 

VIII. DIFFUSION LAW AND ACCOMPANYING GAS STATISTICS 

Combining the generating virial relation (f5"Tj) with (1431) and (|4"5)l , we can write 

(exp {J n(t, x)ln[l + a + jj(x)} dx } ) 

F (t,R,V;n) — t — = (52) 

exp (n J a g[x) dx ) 

l / R,V, (1+ct) n 



= F (t, i?, V; (1 + er) n) <^ exp j y x) In [1 + ip{x)/(l + a)} dx 
As the consequence, after special choice ip(x) = 4>(p) and integration over V , we have 



(exp{/n(t,p) ln[l + a + ^(p)]d 3 p}) fl n 

W(i,i?;n) r f ^rr 3 t — = (53) 

exp \nj a F{p) d 6 p } 



W{t, R; (1 + a) n) { exp <^ / n(t, p) In [1 + 0(p)/(l + tr)] cf/o 



fl, (l+<r)n 
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This generating DVR directly connects the diffusion law W(t,R;n) = Wo(t, R;n) — J Fo(t, R, V; n) d 3 V with gas 
statistics in configurational space as described by one more GF 



= W +jr ^ /•••/ W s <f>{ Pl ) . . . <f>(p s ) 
s=i s - Jl Js 

= W (t, R- n) ( exp j / n(t, p) In [1 + <j){p)] d 3 p 



(54) 



R,n 



where now J\ . . . = J . . . d 3 pj . Clearly, the integrals in both numerator and denominator on the left in (|52|) and (|53|) 
are formally diverging at a ^ , but the divergencies definitely compensate one another, due to equalities (l4l?]) - (l50"l) . 
since at p — > oo random fields n(t, x) and n{t, p) behave like equilibrium Poissonian ones. 

Because of arbitrariness of a and ip( x ) or 4>(p) > each of relations (|52|) and (f53| produces two independent relations. 
Namely, for example, ([53]) yields 



W(t, R; (1 + a) n) 
W(t,R; n) 

or, writing a = exp (£) — 1 , 



exp 



[ n(t, p) In (1 + a) - n E(p) a } d 3 p 



R, n 



W{t,R- e^n) _ ( e*P { U n{t, p) d 3 p } ) R n 
W(t,R;n) ~ (e^{an(t,p)d 3 p}) e n q 

And, taking <j>(p) = (1 + a) [exp (tp(p)) - 1] = e« [exp (cp(p)) - 1] , 
(exp{Jn(t,p)[£ + <p(p)]d 3 p}) Rn 



exp 



n{t,p)(p(p)d p 



{exp{Jn(t,p)£d 3 p}) R n \ k j J / R,ein 

The first of these three equalities is the same as 

W(t, R; (1 + a) n ) = C{t, R, a ; n} , 
which follows from (fBTj) at ip(x) =0 after integration over BP's velocity, with GF 

C{t, R,4>;n] = J C{t, R, V, ip(x) = <j>(p) ; n} d 3 V 

W{t,R, <j>;n} 



(55) 



(56) 



(57) 



(58) 



C (t,R;n) + f^^ J C s (t,R,p^;n) 0( Pl ) . . . <f>(p s ) dp 



exp [n f E{p) 4>{p) d 3 p] 



(59) 



of CFs integrated over all velocities, C s (t,R,p^;n) = f J C s dv^d 3 V. 

Now, let us discuss the above approximate asymptotic expressions for W(t, R;n) at / = in the light of the exact 
relations ([55 ]) - ([55 ) or ([55])) and ([57]). 

In the Gaussian model (|20[). with -D(n) oc 1/n, and in our approximation (l23l) left side of (|55[) and ([56]) looks as 



exp 



(1 + cr) 3/2 exp(-az) 



(e<-l) 



(G) 



Z)(cxp (£) n 



(t,i?) 



W(t, -R; (1 + g) n) 
W(t,R; n) 



{ y ( Av ' ^ ) 



1 



1 + z 
I 



(1 + <J) 3 ' 2 
(l + az/{l + z)) AN -+ 3 / 2 
/ * is]\ W(t,R; e*n) 



(60) 



W(t, R; n) 



(61) 



where z = R 2 /4D(n)t cx n and cr = exp £ — 1 . 



At z = both these expressions reduce to (1+cr) 3 ' 2 
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cxp [(3/2) £] . From viewpoint of right-hand side of (j55|) - 
([56]) it says that the integral J n(t, p) d 3 p = N(t) , - 
representing number of atoms in some relevant volume 
around BP, - consists of non-random constant 3/2 and 
random component which behaves exactly like N(t) in 
equilibrium, with the same Poissonian statistics. Sym- 
bolically, N(t) = 3/2 + N eq (t) . 

Literal interpretation of this statement would mean 
that, strangely, N(t) takes non-integer values. How- 
ever, such interpretation would be ill-advised, since the 
proportionality D(n) oc 1/n, - which just causes fac- 
tor (1 + af/ 2 = exp [(3/2) f] (through W(t, 0; n) cx 
(AitD(n)t)~ 3 '/ 2 ), - is approximate and applicable only 
when 

(n^Vot)- 1 < n < (47ra 3 /3) _1 

and, besides, (1 + a) n = exp (£) n also lies within these 
bounds. What we can state, is that contribution 3/2 
in N(t) — 3/2 + N eq (t) indicates BP-atoms statistical 
correlations produced, of course, by (actual or virtual) 
collisions of atoms with BP. At z = , in that way, that 
are naturally positive correlations: if in spite of t r 
BP has not gone away far from its start position, then 
number of atoms in its surroundings on average is greater 
than under initial equilibrium (for fixed BP), so that 

AiV(i|0,n) = (N(t)) R=0 , n - {N eq ) n = 3/2 (62) 

with {N eq ) n = (N(t))% . 

At finite z > , obviously, in both the cases (pO)) 
and (fFJTj) AN(t\R, n) decreases and becomes negative 
at z > 3/2 and z > 3/2AiV_ , respectively, again in 
agreement with intuition: the greater is current BP's 
distance from its start point, the smaller is number of 
atoms around it, i.e. "number of obstacles" to its flight, 
in comparison with equilibrium. But numerical and sta- 
tistical characteristics of this atoms' "shortage" in cases 
(f60|) and (j6"Tj) are qualitatively different. 

Namely, the Gaussian model, according to (jSUl or 
(f2Tj) . must be associated with infinitely growing short- 
age, AN(t\R,n) = 3/2- z. Thus, if \R\ ~ V Q t then the 
shortage is of order of z ~ t/r, i.e. as large as a num- 
ber of missed BP-atom collisions, as if long ago missed 
atoms somehow determined current shortage. Indeed, 
this is physically rather absurd picture! 

Another picture, which arises under our approxima- 
tion, looks much more likely. Here, according to 
and (|6ip . mean shortage is bounded, 



AN = - 
2 



AN - + 2 



1 + 2 



-AN- (63) 



Besides, it obeys essentially non-Poissonian statistics, 
which means presence of statistical correlations between 
composing atoms. One may describe this statistics by 
treating the "microscopic gas density" n(t,p) as "twice 
stochastic" point field ( "spatial random point process" ) , 



i.e. Poissonian field with fluctuating intensity. Then, 
applying in (|6ip integral expansion 



exp 



AN- + - ) In 



1 



z 



exp 



l + z 

z 



-u - 



l + z 



e«-l) 
(e«-l) 



X -^t-t-t ... ' du , (64) 



we can treat factor 



AJV 



T(A7V_ + 3/2) 



3/2 - uz/(l + z) 



(65) 



(with u obeying the gamma distribution) as fluctuating 
excess or shortage (dependently on A7V 's sign) of num- 
ber of atoms in correlation volume near BP. At that, its 
average value equals to (I63|) while most probable value, 
AN m , p Xz) , is smaller, 



AiV m . p . = - 



AA_ + - 



l + z 



-> 1 - A7V_ (66) 



(since most probable value of u is AA_ + 1/2 ). 

Let us use it for estimation of limits in (1551) and (fBUl) , 
that is characteristic correlation number A7V_ . Intro- 
duce quantity 



AM = M 



, AN 



(67) 



It represents excess (or shortage) of total mass (of BP 
and atoms) located inside the correlation volume fi c , 
in comparison with mean mass mnfl c of other regions 
having such volume but not containing BP. Notice, first, 
that in perfectly equilibrium statistical ensemble, - where 
BP's position is fully uncertain, W(R) — I/O, - the 
same quantity mnfl c gives mean mass of arbitrary re- 
gion with volume fl , even though it may contain BP. 
Second, since in our ensemble the limit value of AiV at 
z — > oo appears insensible to z , it in fact also corre- 
sponds to fully uncertain BP's position ("at infinity"). 
Therefore we can expect that the excess-shortage mass 
AM in our ensemble at z — > oo becomes the same as in 
perfectly equilibrium ensemble, i.e. turns to zero. Then, 
if addressing this requirement to most probable value of 
AM, i.e. identifying AN in (|67p with AN m . p .(oc) from 
(1661) . we come to equality 



M + (1 - AA_)m = 



(68) 



which suggests for A N— value A N- = 1 + M/m . 

This value coincides with what was found formally, un- 
der the collisional approximation, in [3 for BP in 
non-ideal gas. Of course, our derivation of (|68j) is less 
formally grounded, and factual value of AiV_ may differ 
from 1 + M/m. Nevertheless, undoubtedly, this result 
correctly reflects role of mass ratio M/m in construction 
of the law of diffusion, and confirms AN- 's indifference 
to gas density. 
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As the consequence, we can improve estimate of the 
correlation volume f2 c and make estimate of its ana- 
logues for higher-order CFs. With this purpose let us 
return to formula (155)) . Due to relations (fT2)l we can 
write 



k=Q 



Ch 



where C s = C s (t,R;n) are CFs smoothed over small 
enough . Combining this with (J39| and DVR ([13]) . 
we obtain 



r(Aiv_) 
r(Aiv_ + s ) 



(69) 



The factor before integral here just has meaning of in- 
verse 3s -dimensional characteristic volume occupied by 
(s + l)-order irreducible correlation between BP and s 
atoms. Thus the 3s -dimensional volume itself is 



AA_ 



(AM 



8-1) 



(70) 



in particular, f2 c = fl c \ — AjV_/n . 

This expression, being supplemented with above es- 
timate AA_ = 1 + M/m , acquires visual explana- 
tion. Namely, at M/m < 1 we see from (|70|) that 

rfJs ~ ((s + lyy/'/n - (s + l) /en . Hence, 3D correla- 
tion volume of a particular BP-atom link (in pj -space) 
grows proportionally to s , that is correlation between 
BP and one of atoms is mediated and lengthened by all 
s — 1 others. Physically, this means that all they are 
participants of a same connected cluster of s (actual or 
virtual) collisions (or "encounters" [IH, [HJ]), and most 

lis 

long-range contributions to fi c s come from cylinder-like 
regions ("collision cylinders") with cross-section ~ tto 2 
and length ~ sA (thus, with volume ~ sira 2 \ = s/n). 

If M/m ^> 1 , then formula (1701 implies, at not too 

large s , value O^J s s ~ Q c ~ M/mn . Its approximate 
independence on s says about relative weakness of (BP- 
mcdiatcd) inter-atom correlations inside the clusters (if 
BP was immovable, they would vanish at all). 

Physical meaning of the proportionality Q c oc M/m 
at M/m 3> 1 also is quite transparent. Indeed, f2 c is 
formed by "collision cylinders" with length A of order of 
atom's flight path during BP's velocity relaxation time 
r, i.e. A ~ r^/T/m. At that, r ~ MA/VTto [II]. 
This value follows e.g. from phenomenological equation 
for BP's momentum relaxation in light dilute gas: 

MV « -{V/\)VTm~ . Hence, A - (M/m) X , thus 
explaining why Q c — AA_/n m ira 2 K m M/mn . 

Above reasonings demonstrate that characteristic "cor- 
relation volume" can be thought as a weighted sum 
(union) , or statistical sum, of variously oriented (and may 
be coupled) "collision cylinders" . That is why volume f2 c 
is determined by 1/n , but not A 3 ^> 1/n , although its 



greatest liner size is determined by A . Such construc- 
tion of the correlations implies strongly non-uniform dis- 
tribution of "density of correlations" inside fi c . For in- 
stance, - as was suggested already in [2?J (and discussed 
also in 5]), - schematically, C\(t, R, p;n)/C\(t, R;n) ~ 
(a 2 /4|p| 2 ) exp(— |/9|/A) at \p\ > a, which reveals hid- 
den motive of relations (l6"9l - (|70l) . At this point, however, 
we enter to "terra incognita" . 



IX. CONCLUSION 

We have tested one more novel and unprejudiced ap- 
proach to the problem about molecular Brownian par- 
ticle (BP) in ideal gas, basing on explicit density ex- 
pansions of (time-depending non-equilibrium) probabil- 
ity distribution functions and on exact dynamical virial 
relations (DVR) obtained directly from these expan- 
sions. Previous approaches were based on the "gener- 
alized fluctuation-dissipation relations" (FDR) [(| 0, [27j 
(about FDR themselves and their other applications see 
e.g. [l3| and references therein), on the BBGKY equa- 
tions and Bogolyubov's generating functional equation 
0-0, H, H3] , on the "stochastic representation of deter- 
ministic interactions" [29j] and path integrals [l], [2j| . So 
much attention to seemingly too particular problem was 
caused by understanding that in fact it is very non-trivial 
for real "honest" statistical mechanics (SM). Though it 
is quite trivial for "Boltzmannian version" of SM ex- 
ploiting "Bernoullian way of thinking" , i.e. belief that 
many-particle chaos can be divided into "independent" 
elementary events with strictly certain probabilities. In 
the framework of this fantasy, long-range statistics of 
random walk of molecular-size BP does not differ from 
standard mathematical Brownian motion (diffusive pro- 
cess, Wiener process, etc.), while in real SM one reveals 
random walk without certain diffusivity and mobility or, 
in other words, with 1/f-type fluctuations of diffusivity 
and mobility. This difference, firstly, illustrates Krylov's 
statement [25[ that the only "elementary events" what 
generally can be independent and have certain proba- 
bilities are whole phase trajectories of one or another 
many-particle system. Secondly, highlights where one 
should search for origin of various 1/f noises. They, along 
with accompanying historical correlations 0, manifest 
uniqueness of any factual phase trajectory (experiment). 
Thirdly, the mentioned difference once again shows that 
"surprises in theoretical physics" do continue, and failure 
of the "law of conservation of slopiness" 30] (in molec- 
ular chaos considerations) supplements their collection. 
Hence, the problem under our attention is of principal 
importance. 

We developed analysis of our present system in terms of 
conditional statistical characteristics of microscopic gas 
density under given BP's path passed during total ob- 
servation time. We transformed the DVR into simple 
and intelligible relations between these gas characteris- 
tics and BP's diffusion law. First of all, principal re- 
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lation between its scaling exponents in respect to mean 
gas density and BP's path value, on one hand, and con- 
ditional mean number of gas atoms covered by instant 
statistical correlations with BP, on the other hand. We 
demonstrated that such relations predict essentially non- 
Gaussian BP's law of diffusion possessing nearly power- 
law long tail. Moreover, they, - in combination with nat- 
ural heuristic reasonings, - allow quantitative estimate 
of the tail's exponent even without formal summation of 



the explicit expansions. 

The ancient classical "kinetic theory of gases" is unfa- 
miliar with so remarkable relations, but now gets chance 
to enhance its "mental outfit" and language and include 
notion of molecular chaos which produces 1 /f noise. Si- 
multaneously, of course, we need in more complete math- 
ematical investigation (if not exact solution) of the prob- 
lem about BP in ideal gas, along with other problems 
about spatial-temporal statistics of relaxation and trans- 
port processes in many-particle systems. 
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